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An 5(7(2) slave boson formulation of bond-type spin nematic orders is developed in frustrated 
ferromagnets, where the spin nematic states are described as the resonating spin-triplet valence 
bond (RVB) states. The d-vectors of spin-triplet pairing ansatzes play the role of the directors in 
the bond-type spin quadrupolar states. The low-energy excitations around such spin-triplet RVB 
ansatzes generally comprise the (potentially massless) gauge bosons, massless Goldstone bosons, 
and spinon individual excitations. Extending the projective symmetry group argument to the spin- 
triplet ansatzes, we show how to identify the number of massless gauge bosons efficiently. Applying 
this formulation, we next (i) enumerate possible mean field solutions for the S = | ferromagnetic 
J1-J2 Heisenberg model on the square lattice, with ferromagnetic nearest neighbor Ji and competing 
antiferromagnetic next-nearest neighbor J2, and (ii) argue their stability against small gauge fluc- 
tuations. As a result, two stable spin-triplet RVB ansatzes are found in the intermediate coupling 
regime around Ji : J2 ~ 1 : 0.4. One is the Z2 Balian-Werthamer (BW) state stabilized by the Higgs 
mechanism and the other is the 5(7(2) chiral p-wave (Anderson-Brinkman-Morel) state stabilized 
by the Chern-Simon mechanism. The former Z2 BW state in fact shows the same bond-type spin 
quadrupolar order as found in the previous exact diagonalization study [N. Shannon et ah, Phys. 
Rev. Lett. 96, 027213 (2006)]. 

PACS numbers: 



I. INTRODUCTION 

Recent theoretical progress has revealed that a cer- 
tain class of frustrated nmgnetsi^^£& 7 -&2i u i shows 
spin nematic state a 1 ' 11 as their magnetic ground states, 
where the spin quadratic tensor, Kj\^ v = (Sj^Si^) — 

^(Sj ■ Si) with {n,v — 1,2,3), exhibits a long-range 
order, while the spin moment (Sj^) remains disordered. 
Such spin nematic states can be classified into the chi- 
ral type (p-nematic) and non-chiral type (n-nematic) 
states^ according to the parity of the spin quadratic 
tensor. Namely, the antisymmetric quadratic tensor 
Pjl,X = ex^Kji^u is nothing but the the vector chi- 
rality, while the symmetric part — non-chiral one - 
plays the role of the spin quadrupolar moment, Qji^ v = 
^(Kji^^ + Kji^^). The latter ordered state is a spin ana- 
logue of the nematic state well known in liquid crystals^ 
where the order parameter is characterized by the so- 
called 'director vector' d(r) in the form 

Qfiu(r) = d„(r)d„(r) - h^\d{r)\ 2 . (1) 

From this analogy, the spin quadrupolar states are often 
dubbed simply as the 'spin-nematic' states^ 

Depending on how the spin quadrupolar mo- 
ments are microscopically organized, spin nematic 
states have two distinct classes; (i) site-type ne- 
matic statesii±^i±£i±k±L±&±£ and (ii) bond-type ne- 
matic states ±££d£&2&iiL12. The former types of nematic 
orders are realized in the spin 1 bilinear-biquadratic 
model, H s =i = J2 {ij )[JSi ■ Sj + K{S l ■ Sj) 2 }, where 
the quadrupolar moments constituted at respective sites 
exhibit the long-range order due to the strong bi- 
quadratic coupling ] 11 i 12 i 13 i 15 i 17 Ground state wavefunc- 
tions of these site-type nematic states can be es- 



sentially factorized into decoupled 'vacuums', which 
arc defined on respective sites. Thus, their spin- 
wave theorie s 12 i 13 ' 15 i 16 i 18 i 19 including low-energy effective 
theories 2 ^ were well-established. Namely, the elementary 
excitation around such a site-factorized vacuum is also 
given by a linear combination of bosons introduced at 
respective sites. 

The simplest localized spin models which allow 
the second class of spin nematic states - bond-type 
nematic states - are the spin one half frustrated 
ferromagnets fLL2d£&iZ&&12. which could be realized in 
a certain family of layered cuprate o 21 i 22 ' 23 i 24 ' 25 and 
vanadate o 26 ' 27 and also in solid 3 He films.— For example, 
in (CuX)LaNb 2 7 (X=Cl,Br),2Ii22 Cu 2 + ions, having a 
localized spin i, compose a square lattice, while the anion 
X~ locates at the center of the square instead of the bond 
center. As a result, the nearest neighbor (NN) exchange 
interaction J\ between the localized spins becomes ferro- 
magnetic because of the Goodenough-Kanamori rule, 29 
while the next nearest neighbor (NNN) interaction Ji be- 
comes antiferromagnetic; the modcl-Hamiltonian is given 
by 

H = -JiJ2 Sj • S ; + J 2 ]T ■ S h (2) 

(3,1) am 

with Ji, J2 > 0. The preceding exact diagonalization 
(ED) studies for this spin one half square lattice J\- 
J2 model 5 indicated that the d-wave bond-type spin ne- 
matic order develops in the intermediate parameter re- 
gion, Ji ~ 2J2. Namely, strong ferromagnetic exchange 
interactions favor the spin-triplet valence bond forma- 
tions between two neighboring spin one halves, while, 
simultaneously, these two spin one halves try to change 
their partners quantum-mechanically by way of the NNN 
antiferromagnetic exchange interactions. This leads to a 
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kind of resonating spin-triplet valence bond state, where 
the quadrupolar moment organized at each neighbor bond 
exhibits the following antiferro-type configuration with 
the uniform amplitude; 

Qo - ,j+a>,22 " Q(j,j+x),n = Q{j,j+f»,n - Q(j,j+$),22 > 0(3) 

Similar bond nematic order phases were also found in 
other frustrated ferromagnets, such as a zigzag spin 
chain 2-6,7 containing ferromagnetic J\ and a triangular 
lattice multiple-spin exchange modeli&SJii 

In contrast to the site-type nematic states, however, 
when attempting to construct a mean-field description 
of these bond nematic states [as well as their spin wave 
theories], one could immediately reach a more funda- 
mental question; how their ground state wavefunctions 
themselves should be described? Namely, since a sin- 
gle spin one half at each site is supposed to participate 
equally in the spin-triplet formations on its four ferro- 
magnetic bonds [in the square lattice case] , their ground 
state wavefunctions are no longer described by any kind 
of 'site-factorized wavefunctions'. 

In this paper, we will construct an SU(2) slave-boson 
mean- field theory of the bond- type spin nematic states, 
which are described as the resonating valence bond 
(RVB) states of the spin-triplet bonds. After splitting the 
original spin operator into the bilinear of the spinon fields 
(fermions)^^^ Sjli = ± /]>,,] we first in- 

troduce the spin-triplet pairing ansatzes into the ferro- 
magnetic exchange bonds as 

Ev,M = (fL[*Mil>h ( 4 ) 

Aj.ai = (fiaiio-20-vUpfjp), (5) 

where Dy (Ejj ) describes the d- vector of the spin-triplet 
pair condensation 36 ('spin-orbit' hopping integral) . In 
fact, these two-types of the d-vectors, i.e. that in the 
particle-hole channel and in the particle-particle chan- 
nel, precisely mimic the director vector d(r) of nematic 
states in liquid crystals [see Eq. |(TJ)]; in the mean- field 
approximation, the quadrupolar order parameter is given 
by 

Qjl,nv — — ~^{Eji,iiE* lv — —5^ y \Eji\ 2 ^ + h.c. 

-\{D 3 L,D* lu - l -5, v \D 3l \ 2 ) + h.c. (6) 

Moreover, the vector chiral order parameter is given by 
the products between these two <i-vectors and their re- 
spective spin-singlet ansatzes in the form 37 

™ji,A =0 W'-^ji.A ~ Xji E ji,\) 

1 

~ 2 Vl3i D *ji,\ ~ VjiDji,x) , (7) 

where Xjl iVjl) stands for the spinless hopping integral 
(spin-singlet pair condensation) ; 30 i 31 1 33 

Xjl = {fjafj«), Vji = (fj a [(-i)o- 2 ] a pfit3)- (8) 



Thus, one can naturally employ the spin-triplet slave bo- 
son theory as a mean- field description of the spin nematic 
orders. 

In Section II, we will introduce an SU (2)-formulation 
of the spin-triplet mean-field ansatzes, where we exten- 
sively use the 2x2 matrix representation originally intro- 
duced by Affleck et aL,— instead of the usual Nambu vec- 
tor. This representation [see Eqs. (fT0|) and (fT2|) ] clearly 
dictates that the low-energy excitation around any spin- 
triplet RVB state generally consists of (gapless) Gold- 
stone boson and (potentially gapless) gauge boson. It is 
widely known that the existence of the gapless gauge fluc- 
tuations is crucial to the instability of the starting mean- 
field ansatzesi 34 i 35 ' 38 Thus, we will next argue the spin- 
triplet extension of the projective symmetry group (PSG) 
arguments. Without resorting to any microscopic calcu- 
lations, this extension enables us to identify the number 
of the massless gauge bosons for any given mixed ansatz 
having both spin-triplet and spin-singlet link variables. 

Armed with these general formulations, we study in 
Section III the ferromagnetic J\-J% Heisenberg square 
lattice model defined in Eq. @ , thereby finding two sta- 
ble spin-triplet RVB ansatzes in the intermediate cou- 
pling region, J\ ~ 2J 2 . One is the Balian-Werthamer 
(BW) type triplet pairing stated having the coplanar 
configurations of the d- vector, d(k) oc xk x + yk y , while 
the other is the chiral p-wave stated having its d-vector 
all pointing in the same direction d(k) oc z{k x -(- iky^ [see 
Fig.[Hb)]. The PSG arguments indicate that, in general, 
all the non-magnetic (gauge) excitations in the BW state 
have finite Higgs mass. Thus, this ansatz — Z 2 BW state 
— is stable against any type of small gauge fluctuations. 
On the other hand, the chiral p-wave state does not break 
any of the SU{2) gauge symmetry. Instead, it breaks the 
time-reversal symmetry and all the mirror symmetries. 
As a result, nonmagnetic (gauge) bosons are endowed 
by the Chern-Simon term with the topologically-induced 
mass. Thus, this SU(2) chiral p-wave state is also stable 
against any small gauge fluctuation. Though both the 
BW and chiral p-wave states exhibit spin quadrupolar 
orders, the BW state especially shows the same config- 
uration of quadrupolar moments as the bond-type spin- 
nematic order found in Ref. la. Hence, we further discuss 
possible experimental features of this BW state, mainly 
focusing on its magnetic excitations. 

Section IV is devoted to the summary and open issues. 
The relation between our Z 2 BW state and the time- 
reversal topological insulator recently discussed in the 
various literature o 41 i 42 ' 43 i 44 ' 45 is briefly mentioned. We 
also propose those combinations of the triplet and singlet 
ansatzes which describe the vector chiral order having 
no finite director vector,— i.e. Pji.ji ^ and Qji,^ = 0. 
Those readers who want to make the SU (2) slave-boson 
study in frustrated ferromagnets to be a controlled anal- 
ysis might as well consult the appendix A, where we de- 
scribe the large- AT generalization of frustrated ferromag- 
netic spin models. 
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II. S(/(2)-FORMULATION OF SPIN-TRIPLET 
RVB STATE 

A. Matrix representation 

The slave-boson formulation begins with describing the 
spin operator by the bilinear of fcrmion fields; 2Sj^ = 
fj a { a ^] a PfjP- The en l ar g e( i (fermion's) Hilbert space re- 
duces to the physical (spin's) Hilbert space, provided that 
the following local constraints are strictly observed at 
each site: 

/jaMa/3/,73 = 1, ftffjl = /ji/j'T = °- 

In the partition function, these local constraints are im- 
plemented as the coupling between the fermion (spinon) 
fields and the temporal SU(2) gauge fields a" r {v — 
1,2,3)^31^ 



da T d& d^t exp 



/3 



drC 



(9) 



j v=l 



where and stand for the 2x2 matrices 



/j,T 



fh hi 

"-fat 



(10) 



The spin Hamiltonian part TL becomes quartic in the 
fermion field (\I/-field). Depending on the sign of the 
exchange interaction, we decompose this quartic term 
into the Stratonovich-Hubbard variables in two alterna- 
tive ways; 



dU Bin dU tTi da T d¥d$> exp 







drC 



(11) 



1 



C =2 E ^[*5(^0 + E ia lr^i] 
3 



\ 2 -\n 



|E{(-i^i 2 -^H 2 ) + Tr[*t^ ; ]}. 



only over the d-vectors in the excitonic/Cooper channel. 
In fact, this integration precisely reproduces the ferro- 
magnetic exchange interaction, 



iSj-St = -^2(fj a [cr»] a 0fi0)(fh[<T»] 7 6fjs) 

3 



while that over the singlet variable leads to the antifer- 
romagnetic exchange interaction ) 30 i 31 i 33 i 34 ' 35 

AS, ■ S l = -(flfla)(f! f j0 ) 

Thus, the slave-boson formulation of mixed Heisenberg 
magnets generally requires us to use the spin-triplet 
link-variable UjJ 1 for every ferromagnetic bond and the 
spin-singlet link variable C/J ; m for every antiferromagnetic 
bond. 

The saddle point solutions of Eq. (fT2")) lead to the cou- 
pled gap equations for these link- variables, i.e. Eqs. ([!]), 
(O, and (JSJ), whose right hand sides are self-consistently 
given by these mean-fields themselves. In terms of f/jfV 
and Ujl n thus determined, the spin quadrupolar moment 
and vector chirality are given by 

c 3 

-2Q AM „ = Tr[U^ M U^ u ] - -f J2 nuftxUltxh (14) 



-2iP jltX = Tr[^f Uff x ]. 



(15) 



Comparing Eq. (fT2|) with Eq. (fT5|) . notice that the present 
J\-Ji model can have spin quadrupolar order on ferro- 
magnetic bonds, but cannot have vector chirality on any 
links, since Ufj n Ujf\ = 0. Within our formalism, a naive 
mean-field description of vector chiral orders becomes 
possible only in those spin models having either symmet- 
ric anisotropic exchange interactions or antisymmetric 
anisotropic one. In the next section, without making any 
distinction between the n-nematic states and p-nematic 
ones, we will widely call those mean-field ansatzes having 
both finite triplet ansatz and singlet ansatz as spin-triplet 
RVB states. 



B. Low-energy excitations around spin-triplet RVB 

(12) states 



Namely, the triplet and singlet link- variables, 

D 



U sin 



X ij 



-Xij 



ui 



E*. 
-D* E 



(13) 



are introduced as the auxiliary fields for the ferro- and 
antiferro-magnetic bonds, respectively. This is simply be- 
cause the sign of the ferromagnetic exchange interaction 
generally allows us to perform the gaussian-integration 



To see the low-energy excitations around the spin- 
triplet RVB ansatzes, let us first express the spin op- 
erator in terms of the 2x2 matrix representation^ 
Sj^ = iTr^t^ ja^]. Namely, a spin rotation is de- 
scribed by an SU(2) matrix, say hj, applied from the 
left (right) hand side of 
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while physical quantities are invariant under any local 
SU(2) gauge transformation applied from the right (left) 
hand side of 

{'%-i'T\ ■,!,{! z- r r\^- 

For example, both parts of the spin quadratic tensor, 
Eqs. (|14p and (fl~5| . are invariant under this local SU(2) 
gauge transformation. In regard to these two symmetries, 
any spin-triplet mean-field ansatz is generally accompa- 
nied by two types of low-energy excitations: the magnetic 
ones (Goldstone bosons) 47 and the the non-magnetic ones 

(gauge bosons) .^31^35^6,48, 

The former excitations are semiclassically described by 
the deformations of the ri-vectors around its mean-field 
configuration, 

U^^Uf^U±^,r) (16) 

for /x = 1,2,3 with the 3x3 rotational matrix R(x,t). 
Such deformations cost infinitesimally small energy in 
spin models with spin continuous symmetry, provided 
that the variation of the rotation is sufficiently slow in 
space and time. This type of deformations describe the 
Goldstone modes accompanying the spontaneous symme- 
try breaking. 

In addition to this conventional excitation, a certain 
non-magnetic (gauge) excitations also become massless, 
when our starting mean-field ansatz is invariant under 
a continuous gauge symmetry . 30 ' 31 ' 34 ' 35 For example, as- 
sume that the invariant gauge group (IGG) contains the 
[7(1) gauge symmetry {e l6,J3 \8 e [0, 2n)}. Namely, our 
mean-field ansatz is invariant under any rotation around 
the 3-axis in the gauge space, 

^{U^UT}^ 3 = {Vtf,»Ut} ( 17 ) 

for fx = 1,2,3 and dj r = S^a^ T . Then, we can argue 
that the following non-magnetic deformation also com- 
prises the gapless excitation: 

{U;t,Uf}^{U^,U]?}e^, (18) 
a i,r = «|t + a oU,r), (19) 

where a,ji relates to the spatial components of "gauge 
fluctuations" a a (J, r) (a = 1, • • • , d) in the form 

a ji( T ) = U - l) a a a (j,T). (20) 

Specifically, one can expand the effective action in terms 
of these variations a a {j,r) (a — 0, 1, ••■ , d), assuming 
these fluctuations to be much smaller than their units, 
a a(j,T~) <C 2tt. Up to their quadratic order, the effective 



action generally reads as follows: 

d 

F„ = E H^MQKWM-Q) + ••• . ( 21 ) 

a,/3=0 Q 

fl ^) = i|EE e^-^a a (j, t) (22) 

V iuj n q 

with Q = (q,iuj m ). Then, taking into account the U(l) 
gauge symmetry of the mean-field ansatz, one can specify 
the form of the (d + 1) X (d + 1) matrix M(Q), such that 
the quadratic part in Eq. |2T|) reduces to the U(l) gauge 
invariant form as in Eq. (|23j) . 

To see this, introduce the following local U(l) gauge 
transformation in Eq. (fT2| : 

*,-(r) - e-^M^^(r), 

where dj(r) varies slowly in space and time. Under this 
transformation, all changes in the link variables (fT8|) are 
put into the transformation, aji — > aji +8i — 9j and ao — > 
a + d T 8, due to the U(l) symmetry in IGG. Thus the 
effective action around Q ~ is literally transformed as 

F — > 

r gauge ^ 

M ^{Q)(a a + d a 6)(Q)(a + d p 9){-Q). 

a,P Q~0 

However, the free energy should have been invariant un- 
der any gauge transformation, since gauge degrees of free- 
dom can be absorbed into the integral variables, VP-fields. 
This requires that M(Q) must precisely reduce to zero 
at Q — 0, so that the quadratic part of the action takes 
the t/(l) gauge invariant form, e.g. 

^ = i E E ^UQ)fa(-Q) + • ■ • , (23) 

where / 7 = £a/3 7 9 Q a^ stands for the field 
strength i 30 i 31 i 34 i 35 It is well-known that this maxwell 
form does not suppress the gauge fluctuation effi- 
ciently. Especially, when the mean-field ansatz have 
its fermionic excitations fully gapped and when d = 2, 
these massless gauge fluctuations destroy the mean- 
field ansatz itself , 34 ' 35 ! 38 apart from some exceptional 
cases . 49^^,^,53 ,54 F u owm g the literature,— we call 

in this paper such spin-triplet mean-field ansatz as the 
gapped U(l) [or SU(2)] state. 

On the other hand, if the starting mean-field ansatz 
has no continuous invariant gauge group (IGG) like in 
Eq. (|17[) , the local minimum condition imposed on mean- 
field ansatzes generally requires all the eigenvalues of 
M (Q) to be positive. Therefore, all the gauge fields have 
finite Higgs mass around any Q; 

d 

^auge = E E M a (Q)a a (Q)a a (-Q) + ■■■ (24) 

Q a=0 
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with M a (Q) > 0. In contrast to the maxwell form dis- 
cussed above, this finite Higgs mass suppresses any small 
gauge fluctuation completely, so that the starting mean- 
field ansatz is always guaranteed to be (at least locally) 
stable. Such ansatzes are usually dubbed as the Z 2 state. 

The efficient way to confirm the absence of the contin- 
uous IGG was introduced by Wen ; 34 ' 48 where he pointed 
out the sufficient condition for its absence. We can ex- 
tend his argument to the spin-triplet RVB states also. To 
see this, let us begin with the calculation of the SU(2) 
flux defined on a plaquette by multiplying link-variables 
along the closed loop in a regular sequence, where ei- 
ther Uf™ or one of U^ 1 should be chosen on each link. 
For example, when the loop is given by a triangular path 
i — > j — ► k — > i, one can have an SU(2) flux by UijUjkUki, 
which always transforms in a gauge-covariant way; 

UijUjkUki — ► g\ ■ UijUjkUki ■ gi 

under ^ — > Sj^j- As such, the relative angle sub- 
tended by two distinct St/ (2) fluxes derived from the 
same base-site, such as UijUjkUki and UijUjiUu, contains 
non-trivial gauge-independent information, provided that 
the two triangular paths, (ijk(i)) and (ijl(i)}, are differ- 
ent with each other. Note that, even out of the same tri- 
angular loop, we can have two distinct fluxes, when one 
of its three links has two different types of spin-triplet 
ansatzes, Ufj\ ^ ^IjV in this case, we should regard 
that U^UjlU^ and U\f^U'- k U' ki are two distinct fluxes 
obtained from the same base-site i. 

Having all SU (2) fluxes thus obtained in hand, one can 
readily see that, (i) if two distinct SU(2) fluxes obtained 
from the same base-site are not collinear with each other, 
there is no continuous IGG in that mean-field ansatz. (ii) 
If all the distinct fluxes obtained from the same base-site 
are pointed along one direction in the gauge space, say 
along the 3-axis, the ansatz could have a certain U(l) 
gauge symmetry around this 3-axis, just like in Eq. (fT7l) . 
One can also confirm that, (iii) the ansatz can be in- 
variant under a certain SU (2) gauge symmetry [so-called 
SU(2) state], if all the SU(2) fluxes are proportional to 
the unit matrix. 

This 'non-collinearity' argument of the SU(2) fluxes 
concludes the (local) stability of each ansatz against 
gauge fluctuations very efficiently, without resorting to 
any microscopic calculation. Thus, it substantially helps 
us to find a better spin-triplet mean-field ansatz as in the 
case of spin-single RVB ansatzes . 34 ' 48 

III. Ji-J 2 FRUSTRATED FERROMAGNETIC 
SQUARE LATTICE HEISENBERG MODEL 

In this section, we will apply the spin-triplet slave- 
boson mean-field formulation onto the spin-i J\-J 2 
mixed Heisenberg model ([2]) on the square lattice with 
ferromagnetic nearest neighbor (NN) J\ and antiferro- 
magnetic next nearest neighbor (NNN) J 2 . As was de- 
scribed in the previous section, we always decompose the 



ferromagnetic NN bond into the spin-triplet ansatz and 
the antiferromagnetic NNN bond into the spin-singlet 
ansatz. 

A. Mean-field solutions 

To be specific, we have numerically studied the var- 
ious local 'stable' minima of the mean-field free energy 
given in Eq. (|I2[) . assuming that the magnetic unit cells 
(MUC) are either (i) original square-lattice unit cell or 
(ii) 2 x 2 of the original unit cell. The dimension of the 
(real-valued) parameter space in each case becomes (i) 
32(+3) and (ii) 128(+12). Starting from a randomly cho- 
sen initial point in these multiple dimensional parameter 
spaces, we perform the Newton-Raphson method, only 
to reach a certain local minimum of the mean-field free 
energy E mi (per the magnetic unit cell); 

<j/)eMUC 

+ J i E (M 2 + fel 2 ) 

(OO)eMuc 

16tH ^ J J MBZ 

with (i) v = 4 or (ii) v = 16. Here, the summation 
over jl is taken within each magnetic unit cell and A Q 
denotes the spinon energy band. We have repeated this 
procedure from 50 times to 300 times for each parameter 
point, i.e. (Ji,J 2 ) = (sin0, cos#) with < 6 < 5 . In this 
way, we enumerated various spin-triplet RVB ansatzes. 

Throughout this extensive search, we found basically 
three distinct RVB ansatzes having both spin-triplet link- 
variable on each NN bond and spin-singlet link-variable 
on each NNN bond. All of these three do not break any 
translational symmetries of the original unit cell, i.e. T x 
and T v . 



1. Z2 Balian-Werthamer state 

The first one is a sort of the Balian-Werthamer (BW) 
stated where the d-vector on the NN x-link is perpen- 
dicular to that on the y-lmk, 

U {jlj+x) ,„ = i^iDa 2 , U^l j+0) !fi = i5^Da 2 , 
Uf* ]+±±§) = xo-3 ± Wi, ia p = 0. (26) 

'£>', 'x' and '77' above correspond to the real parts of 
Eqs. |5|) and (|HJ), respectively. This RVB state exhibits 
the same antiferro-type configuration of quadrupolar mo- 
ments as the bond nematic state found in Ref.[5j. Namely, 
the nematic order parameters on NN bonds show 

2 I 

Qjl,u = ' Qjl,2i = Qjl,33 = g-D (27) 
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for the cc-direction and 
2 2 

Qji,22 = i 



— Qji,33 — t;D 2 



(28) 



for the y-direction, where Qji^v = for /i ^= v (see 
Fig. [TJ. While this mean-field ansatz breaks the mir- 
ror symmetry P xy which interchanges x-link and y-lmk, 
it is invariant under the following combined symmetry 
and gauge transformations: G X T X , G y T y , Gp^Px, Gp y P y , 
Gpi P' and GtT. The respective gauge transforma- 
tions read 

G x = G y = a , G Px = itrti-lY', G Py = ^(-1)*-, 
Gp iv =ia 2 (-l) jv , G T = (-lf x+iy . (29) 

Here T refers to the time-reversal symmetry, while P' xy 
stands for the mirror symmetry P xy accompanied by an 
appropriate spin-rotation about the 3-axis by 7f/2. 




FIG. 1: (color online) (a) 2z 2 — x 2 — y 2 type quadrupole mo- 
ment formed by two S = 1/2 spins on each bond, (b) J1-J2 
model and the configuration of the quadrupole moments on 
bonds in the Z 2 BW state [see Eqs. <(27j) and ([28)) ]. 



Provided rjx 7^ 0, the ansatz supports two non- 
collinear SU{2) gauge fluxes, 

^'j+£>4^<J+Sj+*+S),2^(T+i+yj> x X°3 + Wi, (30) 

U% j+x) ,^ xJ+i _ y) ^ x _ yJ) oc X a 3 - Wl . (31) 

Hence it is protected from any small gauge fluctuation 
by finite Higgs mass. We call this ansatz as the Z2 BW 
state. The spinon's band dispersion A Q of this Z2 state is 
comprised of two doubly degenerate bands, both of which 
are always separated by a finite energy gap in the entire 
Brillouin zone, [— tt,tt] x [— tt,tt]; 

Ai, 2 = -A 3 ,4 = {A 2 (s 2 x + 4) + B 2 c 2 x c 2 y + C 2 S 2 x s 2 y }H^) 

with (s„,c„) = (sinfcu, cosfc^) and (2A,B 7 C) = 
(JiD, J 2 X, J2V)- 



2. SU (2) chiral p-wave state 

The second ansatz we found is the chiral p-wave 
[Anderson-Brinkman-Morel (ABM)] stated in which all 
the d-vectors on the NN-bonds are collinear, while the 
tZ-vector on the x-link acquires extra phase i in relative 
to that on the y-link, 



7-rtri 



iS ll 3Da2, 



Trtri 



(33) 



Namely, two 'D' appearing in the first line stand for the 
real and imaginary part of the d-vector respectively. Be- 
cause of this relative phase factor, this ansatz has its 
fcrmionic band-dispersion fully gapped in the whole mo- 
mentum space; 

Ai, 2 = ™A 3 , 4 = \k = {A 2 {s 2 x + s 2 y ) + B 2 c 2 x c 2 y y. (34) 



In this state, all NN bonds have the same ferro-nematic 
order Q A33 = -§£> 2 , Qji,u = Qjl,22 = \D 2 . 

The IGG of this chiral p-wave state contains the fol- 
lowing three continuous gauge symmetries: 



{e i{-iy^e<r^ e i(-iy-ea^ e i{-iyve^ ee ^^n)}. (35) 

Correspondingly, the low-energy effective theory in the 
gauge (non-magnetic) part consists of three maxwell 
forms around q = (tt,tt), (tt,0) and (0, n) respectively. 
Namely, above continuous gauge symmetries require that 
the following three types of non-magnetic deformations 
constitute the U(l) gauge invariant effective actions: 



al T = (-iy* +j «ao(j,T), (36) 
a} jT = (- lj^-oo (7, r), (37) 



"3,T 



(-l) J *ao(j,r). 



(38) 



Though these three types of gauge fluctuations are 
not suppressed by finite Higgs mass, the ansatz it- 
self is still protected by the so-called Chern-Simon 
mechanism^5^££i£^ 

To see this, notice that the ansatz (1331) breaks all the 
mirror symmetries P x , P y , P xy and the time-reversal 
symmetry T. Instead, it is invariant only under these 
mirror symmetries accompanied by the time-reversal 
symmetry Gpq- ■ PT or under the spatial inversion sym- 
metry Gp^R^. The respective gauge transformations are 
given by 

Gp v t = to, G n ^ = G PxT = (-iy* +j v, 
Gp IS T^(-a 3 ) J - +3 ». (39) 

This magnetic point group clearly allows the spontaneous 
Hall conductance of the 'spinon', like in the chiral spin 
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statej 51 i 52 ' 53 In fact, corresponding to the three contin- 
uous gauge symmetries given in Eq. (|55|) . we have three 
conserved 'charges', all of which are accompanied by fi- 
nite quantized transverse conductance a xy = w-- As a 
result, the effective actions around q — (n, tt), (0,7r) and 
(tt, 0) acquire the Chern-Simon term in addition to the 
maxwell form^££L£2£L 

-Fgaugc = j ' dx 2 dT-^-a^d v a\e^ v \ + (maxwell form). 

This Chern-Simon term endows the apparently massless 
gauge boson with a finite energy gap^ 



3. Z-2 collinear state 

The third stable ansatz we found is the 'collinear' state, 
where all <i-vectors are pointing to the same direction, 



rrsin 



(40) 



showing ferro-nematic order = — ^D 2 , Qji : ii 

Qji.22 — \D 2 . Although having the same spin- 
quadrupolar moment as the previous one, this collinear 
ansatz is a distinct quantum order state from the SU (2) 
chiral p-wave state. It preserves mirror symmetries as 
well as the time-reversal symmetry. In fact, one can see 
that all the discrete symmetries of the original square 
lattice are recovered, when combined with the following 
gauge transformations: 

G, = G y = a , G Pm = iai(-l) j *, G Py = ^a 1 (-l) J ^ 
G Pxv = l, G r = (-l)^ +i ". (41) 

Having the non-collinear SU(2) gauge fluxes as in 
Eqs. (|30[) and (|31[) . all the gauge fluctuations around this 
ansatz are suppressed by finite Higgs mass. We hence 
call this state as Z 2 collinear state. 



B. Phase diagram 

The mean-field energy for these three ansatzes are 
plotted in Fig. HJa) with (Ji,J 2 ) = J(sin(9,cos#). Let 
us begin with the lowest energy mean-field solution 
in the well-studied limit, J 2 ^S> J\. In the strong 
J 2 limit, our model reduces to the two decoupled an- 
tiferromagnetic square lattice, so that the knowledges 
of the saddle-point solutions in this limit have been 
well-established3^^4^^£6£L££ Namely, the tt-Aux 
state defined on each square lattice, 

^"+^=^±^1, E/ S 1 J+o>,/.=Kr = ( 42 ) 

with \ = ?7, becomes global minimum, when the mag- 
netic unit cell (MUC) is restricted to the original square 



lattice unit cell. On the other hand, when the MUC 
is enlarged up to the 2x2, the global minimum state 
becomes the staggered dimer state introduced on each 
decoupled square lattice, e.g. 

TT sin — TT sin — ->/,r 77 tri — ir,^ — 

u (j-j+£+y) ~ u U+£,j+v) ~ X CT3 ' u (j,j+»),n ~ la 3,r - u 

rrsin _ rrsin _ n 

U (j+x ,j+2x+y) — u {]+2x,]+x+y) u ' 



TJSin _ rrsin _ n 

U U+V,j+x+2y) — u {j+x+y,j+2y) ~ u ' 

rrsin _ rrsin _ n 

U U+x+y,j+2x+2y) ~ u {]+2x+y ,j+x+2y) - u ' 



(43) 



However, using the variational Monte Carlo (VMC) cal- 
culations, Gros and his co-workers^ have demonstrated 
that, when projected onto the original (spin) Hilbert 
space, the 7r-flux state eventually wins over this isolated 
dimer state. In fact, it is well-established^ that the pro- 
jected 7r-flux state gives the second best variational en- 
ergy in the strong J 2 limit (the best variational estimate 
is obtained from the Neel order stated). 

When increasing the NN ferromagnetic interaction J\ , 
a finite spin-triplet ansatz continuously develops on the 
top of this 7r-flux state, while simultaneously the parame- 
ters r) start to deviate from x, i.e., f] ^ X- This leads to ei- 
ther Z 2 BW state or Z 2 collinear state for Q c \ = 0.66 < 9. 
Thus, the transitions from the 7r-flux state to these two 
Z 2 states are both the second order at the mean-field 
level. Energetically speaking, the Z 2 BW state gives 
a slightly lower mean-field energy than that of the Z 2 
collinear state. 

Notice also that these two Z 2 states are clearly pre- 
empted by the staggered dimer state, Eq. (|43|) . at the 
mean- field level [see Fig. Ufa)]. Observing the situa- 
tion in the strong J 2 limit, however, one can naturally 
expect that, when projected onto the physical (spin) 
Hilbert space, both Z 2 states would win over this iso- 
lated dimer state in the case of a finite J\ . Namely, since 
our Z 2 states are constructed based on the decoupled 
7r-flux states [compare Eqs. (|26l40p with Eq. (gl])], they 
would certainly acquire substantial resonance energies in 
the same way as the 7r-flux state does. On the other 
hand, being factorisable, any isolated dimer state cannot 
gain such resonance energies, irrespective of finite fer- 
romagnetic exchange interactions. Moreover, Fig. [2a) 
indicates that the Z 2 BW asatz is quite energetically 
tunable in the presence of the ferromagnetic exchange 
interaction. Thus, we presume that the Z 2 BW state 
finally dominates in this intermediate coupling region, 
9 cl = 0.66 < 6. 

When 9 c2 ee 0.76 < 0, this Z 2 BW state reduces to 
the U(l) state having no finite rj. Namely, with r\ = 
0, two SU{2) gauge fluxes given in Eqs. (J3DJ and (|3"T]) 
become collinear with each other. Simultaneously, this 
U(l) BW state becomes energetically degenerate with 
the SU(2) chiral p-wave state. Namely, both of them 
have precisely the same mean-field band dispersions ±Afc 
[compare Eq. {31} with Eq. ([3_2]) having 77 = 0] . 

This U(l) BW state is destroyed by the infinitesimally 
small gauge fluctuation. Namely, in the absence of finite 
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FIG. 2: (color online) (a) Mean-field energies (per site) of the 
various ansatzes in the S = 1/2 square lattice ferromagnetic 
J1-J2 model. Note that (Ji, J2) = J|(sin 9, cos 6), where the 
energy unit is taken to be |J|. The Blue line (labeled as B) is 
for the Balian-Werthamer (BW) state, which is the Z2 state 
for 9 c i < 6 < 9 C 2 and which reduces to the U(l) state for 
9 C 2 < 9. The green line (labeled as C) is for the Z2 collinear 
state, while the red line (labeled as D) stands for the SU(2) 
chiral p-wave state. The red dotted line (labeled as A) is the 
doubled 7r-fiux state, where both the A-sublattice and the 
B-sublattice support 7r-flux states respectively [see Eq. (|42[) ], 
The Blue dotted line (labeled as E) is for a set of 'flat-band' 
states (-Bfiat = — §|./|sin#), all of which give the same best 
mean-field energy for 9 c s < 9. The green dotted line (labeled 
as F) is for the staggered dimer state (-E2J' ncr = — ||J| cos#), 
where both the A-sublattice and the B-sublattice support 
staggered dimer states respectively [see Eq. (|43p for its ex- 
ample]. These isolated dimer states are known to be over- 
come energetically by the doubled 7r-flux stated when they 
are projected onto the physical Hilbert space. Since the Z2 
BW state is composed on the top of the 7r-flux state, this 
staggered dimer state is also expected to be overcome by the 
projected Z2 BW state, (b) Expected mean-field phase dia- 
gram in the intermediate coupling region. The transition at 
9 c i is the 2nd order, since the magnetic space group of the Z2 
BW state belongs to that of the 7r-flux state. On the other 
hand, the transition at 9 c s is the 1st order at the mean-field 
level, which one can see directly from the Figure (a). 



77, the non-magnetic deformations denned in Eq. (|19p con- 
stitute the following maxwell form around q — (tt,tt), 

-Fgaugc = J dr J d 2 x{ue 2 + ^Kb 2 } + ■■■ 



where e a (a = 1, 2) and b are defined, from Eqs. (|19I20[) . 

as e a (j,r) = (-l)J'«+iv (d T a a - d a a ) and b(j,r) = 



(— iy x+ iy (d 2 ai — dia 2 )- Since the fermionic excitations 
are fully gapped even without 77 [see Eq. ([52"])], this 
maxwell form is free from any dissipation effect,— e.g. 



d 2 k- 



A 2 B 2 c 2 Jl 



44) 



16tt 2 A| 



> 0. 



Having the time-reversal symmetry [see Eq. (|2"9")) ]. the 
massless gauge fluctuation is not suppressed by the 
Chern-Simon term either^ Consequently, infinitesimally 
small fluctuations of this type of gauge fields lead 
the U(l) BW state into a confining phase having no 
gapped free spinon in its excitation. More specifi- 
cally, those space-time instantons (monopoles) which 
are allowed by the corresponding compact QED action, 
J dr J d 2 x{ue 2 — K cos(e ct ^A Q a / 3)}, proliferate in the 
2 + 1 dimensional spaced lowering a certain magnetic 
symmetries enumerated in Eq. (|29|) . 61 To capture the re- 
sulting magnetic space group of the confining phase, one 
generally need to identify the quantum number carried 
by this monopole creation field . 62 > 63 

For 8 c3 = 0.775 < 6>, these two degenerate ansatzes, - 
5/7(2) chiral p-wave state and U(l) BW state -, are fur- 
ther overcome (energetically) by another ansatz, which 
we dubbed as the 'flat-band' states, 



UJ? = 0, 



uf * 0, 



ia V j T ^ 0. 



(44) 



These 'flat-band' states do not have any finite singlet 
ansatzes anymore and keep on giving the lowest mean- 
field energy (E^ t — — |Ji) for the remaining ferromag- 
netic side, f9 C 3 < 8 < 5. However, these 'flat-band' states 
do not necessarily refer to a specific configuration of the 
spin-triplet ansatzes. Instead, they refer to a group of 
the states all of which give precisely the same mean-field 
energy. For example, these 'flat-band' states include the 
following parameterization of the spin-triplet ansatz: 



r t?" n" rv 

I^Uj+z) ^Uj+z) 

[ E {i,3+y) D UJ+y) D '(j,j+v) ] = P ni 
E', A = a n 2 , E'tfj+f)) = ft n :h 



a ni ■ , 



n, 



n 2 



Tin 



"3 



rii = m 1 ■ m 2 = 0, (45) 



where a 2 + (3 2 = ~ 
trary unit vectors that observe Eq 



and rij and rrij can be arbi- 
Here E k» 



and £)', ,. stand for the real part of Eji t)Ji and Dji 



respectively, while E'- x and D'^ are their respective 
imaginary parts. Thus, only the first one is parity even 



while the others are odd, E'J, „ 



and Dji ^ = —Dij ^. Bearing these in mind, one can eas- 
ily see that this mean-field ansatz always gives the two 
doubly degenerate spinon bands, which are totally flat in 
the entire Brillouin zone, 



Ai, 



-A 



3,4 



Ji 
4 ' 



Because of this feature, all the spin-triplet ansatzes pa- 
rameterized by Eq. (|45[) give the same mean-field energy 
(per site) Eg£ = -|Ji. 
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The emergence of these 'huge' numbers of 'fiat-band' 
states in the strong Ji limit reflects the fact that the 
ground-state order parameter of any Heisenberg ferro- 
magnet (total spin moment) and the corresponding spin 
Hamiltonian are simultaneously diagonalizable. When 
projected onto the physical (spin) Hilbert space, we ex- 
pect that these flat-band states reduce to a fully polarized 
state (ferromagnetic state). 

Observing Fig. [2j please notice that our Z2 BW phase 
appears in larger J2 region in comparison with the previ- 
ous ED studies. Namely, Fig. [5] indicates that its phase 
boundaries are given by Ji : J2 = 1 ■ 1.29 at 9 = 9 C \ 
and Ji : J2 = i : 1.05 at 9 — 9 C 2, while d-wave bond 
nematic order phase was found in 0.4 < J2/J1 ^ 0.6 in 
the previous finite-size studies^ This discrepancy sim- 
ply stems from the so-called 'factor 3' difference, often 
encountered between the Hartree-Fock (HF) spin-singlet 
ansatz and the HF spin-triplet ansatz. If one employed 
a more numerics-oriented formulation^ A appearing in 
Eq. (fl"2j) is replaced by while ^ is replaced by ^jja. 
Consequently, we have J\ : J2 = 1 : 0.43 (9 = 9 c i) and 
J\ : J2 = 1 : 0.36 (9 — 9 C 2), which would be relatively 
comparable with the previous ED result. More quantita- 
tive comparison, however, requires the variational Monte 
Carlo studies based on these spin-triplet ansatzes. 

In summary, we have argued that three spin-triplet 
RVB ansatzes — Z 2 and U(l) BW states and SU(2) chi- 
ral p-wave state — become the lowest mean- field states in 
the intermediate coupling region, J\ ~ 2J2 [see Fig.[2jb)]. 
Among them, both the Z2 BW state and the SU{2) chiral 
p-wave state are stable against any (infinitesimally) small 
gauge fluctuation, while in the U(\) BW state the effect 
of gauge fluctuation is crucial, making spinons confined. 
Using Eq. © , one can easily see that the BW states show 
the d-wave bond- type spin quadrupolar order precisely as 
in Eq. ©. 



C. Magnetic excitations in the BW states 

Here we briefly discuss magnetic excitations in the Z2 
BW state. The 'low-energy' excitation around the Z2 
BW state is composed of three parts; (i) gapped non- 
magnetic excitations (gauge bosons), (ii) gapless mag- 
netic excitations (Goldstone bosons), and (iii) gapped 
fermionic ('J'-field) individual excitations. The gapped 
gauge boson plays only a subdominant role in the spin- 
structure factor, while the latter two contribute signifi- 
cantly to magnetic excitations. Up to the Hartree-Fock 
level, one can easily see that the gapped fermionic ex- 
citation constitutes the continuum spectrum above uj > 
max(Ji|Z?|, J2|x|). When one further takes into account 
the random-phase approximation termS ) 65 i 66 the gapless 
bosonic dispersions emerge below this spinon continuum, 
whose low-energy limit can be described by the matrix- 



formed non-linear a model, 

^magnetic = £ TV [A^iT 1 ^] . (46) 



Namely, the 3x3 matrix R is nothing but the spatio- 
temporally varying rotational matrix of the director vec- 
tor used in Eq. (|16[) . The symmetry argument 6 » dictates 
that the diagonal matrices A M generally take the follow- 
ing form: 



{ A T , A x , Ay } 



r 


CO 




Cl 




Cl 


1 




C2 


) 


C3 


) 


C4 






C2 _ 




C4 




c 3 _ 





•,(47) 



where the director coplanar plane was taken to be the 2-3 
plane. In terms of the semiclassical (gradient) expansion, 
one can directly calculate their respective coupling con- 
stants: 



co = 0, ci 



S y C x 



64tt 2 A3 : 



C2 



C3 



C4 



d 2 k A2{s " +S y ) 



d 2 k 



d 2 k 



647r 2 A| 
JiA 2 s 2 x 



8tt 2 



64^ 2 A3 



64tt 2 X 3 k 



with Ji = A- l \^{{d kx n T ){d ku: n) + 2- 1 A^ 1 9 2 i A fe }, and 

n = 2~ 1 \ k ~ 1 (2As x , 2As y , Bc x c y , Cs x s y ). In addition to 
these massless excitations, we could also have several 
gapped ('optical') magnetic modes, provided that they 
are not damped by the spinon individual excitations i 65 ! 66 
One might also expect a certain characteristic behavior 
of the spectral weight themselves. In fact, Tsunetsugu 
et alr^ and Lauchli et al. 19 demonstrated that the spin- 
structure factor in the site-nematic ordered state exhibits 
the vanishing spectral intensities of the Goldstone modes 
around the r-point. 



IV. SUMMARY AND OPEN ISSUES 

In this paper, we have introduced the spin-triplet slave- 
boson formulation as a mean-field theory for the bond- 
type spin nematic state, which was described as the 
spin-triplet RVB state. Namely, the d-vectors of the 
spin-triplet RVB ansatz constitute the quadrupolar or- 
der, while the combination of the spin-triplet and singlet 
link variables on the same link leads to the vector chiral 
order. 

When applied to the S = | square-lattice frustrated 
ferromagnetic Heisenberg model, our spin-triplet slave- 
boson analysis gives two non-trivial stable spin-triplet 
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RVB ansatzes in the intermediate coupling region around 
Ji : J 2 ~ 1 : 0.4. One is the Z 2 BW state stabilized by 
the Anderson-Higgs mechanism, while the other is the 
SU(2) chiral p-wave state protected by the Chern-Simon 
mechanism. Our slave-boson analysis also found an un- 
stable U(l) BW state as a mean-field solution, which 
possibly gives a route to the realization of spinon con- 
fined quadrupolar ordered states with a certain symme- 
try reduction. The projective symmetry group of the Z 2 
BW state as well as the U(l) BW state is consistent with 
the magnetic space group of the d-wave bond-type spin 
nematic state discussed in Ref. y5|. Both of them exhibit 
the antiferro-type configuration of the bond quadrupolar 
moment shown in Fig. [TJ 

Contrary to a naive expectation, our BW state is clas- 
sified into a 'weak topological (ordinary) insulator' in- 
stead of the 'strong topological insulator' defined in the 
recent literatures! 41 ' 42 ' 43 ' 44 ' 45 Physically speaking, such 
a 'weak topological insulator (WTI)' is accompanied ei- 
ther by no spinon edge states at all or by even numbers of 
the helical edge states. To see that it is indeed a 'WTI', 
one can first deform this ansatz into the U(l) ansatz 
(ji — > 0). Since the fcrmionic dispersion remains gapped, 
the Z2 topological index associated with the filled spinon 
ban d 41 ' 43 is also unchanged. After reaching the simpler 
U(l) ansatz, let us then utilize the Fermi surface argu- 
ment recently introduced by Sato^ His argument relates 
the Z 2 topological index in the superconducting state 
(D ^ 0) with the Fermi surface topology in the corre- 
sponding 'normal' state (D = 0). That is, if a Fermi 
surface in the normal state surrounds odd/even num- 
bers of the time-reversal invariant momentum (TRIM) 
points, the BW state constructed on top of this normal 
state is accompanied by non-trivial/ trivial Z2 topological 
index. Since our normal state is composed of two decou- 
pled u-RVB states at rj = 0, the resulting Fermi surface 
clearly surrounds two time-reversal symmetric ^-points, 
i.e. (0,0) and (ir,n). Thus, our Z 2 BW state should be 
classified into the WTI' (- Z 2 even class -). 

In the followings, we will enumerate several open issues 
and possible extensions of the current work. The most 
immediate open issue is to identify the magnetic space 
group of the confining phase proximate to the 17(1) BW 
state based on the monopole field studies! 62 ' 63 Namely, 
such an analysis gives several complementary informa- 
tions to the direct ED studies of the original spin models 

The fate of the SU(2) chiral p-wave state observed at 
9 C 2 < < 8 C 3 is not so clear either, although we have 
argued its stability against any (infinitcsimally) small 
gauge fluctuation. Namely, previous exact diagonaliza- 
tion studies of the SU{2) spin model did not find any 
T-symmetry breaking ferro-nematic states between the 
d-wave bond-nematic state and ferromagnetic state. In 
fact, it is also possible that, when projected onto the real 
(spin) Hilbert space, the strong gauge fluctuation could 
wipe out this time-reversal breaking ansatz. 

Though we have mainly discussed the quadrupolar or- 
der in this paper, our formulation can also describe vec- 



tor chiral order having no quadrupolar moment^ i.e. 
Pjl,X 7^ and Qji,fiu = 0. In fact, such vector chiral 
order state was observed in the spin one half frustrated 
Heisenberg model having the ring exchange coupling. 4 
When applying the current spin-triplet slave-boson for- 
mulation onto these quantum spin systems, one could use 
the following mean-field parameterization: 

[ E "l D "l D> l ] = [n, n 2 n 3 ], 
[ A -V'jl *fjl ] = [ 7i 72 73 ] , 

where {ni,ri2,ri3} are the normalized unit vectors or- 
thogonal to one another. Namely, such an ansatz gives a 
finite vector chirality, Pji = 2i^2, , 7 q ti q , without any 
quadrupolar moments. We generally have three alterna- 
tive ways to parameterize this vector chiral order, 

l E Jl D "l ^] = ["! ri2 n 3 ], 

[ X", -Vji v"i ] = [ 71 72 73 ] , 

or 

[ E "l E '3l ^] = [«1 ^ »3], 

[ X' jt X"i V"i ] = [ 7i 72 73 ] , 

or 

[ E "l D 'h E 'fl] = [^ ^ 

[ X'ji -V'ji X"i] = [ 7i 72 73 ] ■ 
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APPENDIX A: LARGE- iV FRUSTRATED 
FERROMAGNETIC MODEL 

The mean-field analysis described in this paper be- 
comes exact in the large N limit of the following action; 



element of our Hilbert space: 



N 

{ E f& } |p h y) = o. v i- m = i, 2, a 



o=l 



/3 



Z=y d[/ sm dC/ tn da r ^ aT d* a exp 

£ 4 E Tr t<(^o • 



, (Al) 



The density and spin operators defined in Eq. (|A4|) in 
fact act within this physical Hilbert space. Moreover, 
they observe the following commutation relations: 69 



I 2 -I^| 2 ) 



l^l 



Tr[*f^*^ 



(A2) 



} 



where the summations with respect to the fermion's 
species index a (= 1,- • • , N) were made implicit. The 
integration over the auxiliary fields leads the following 
large- N spin Hamiltonian for frustrated ferromagnets: 



H 



N 



h e {sf ■ s» a + ^v? q } + § E s f ■ s " a - 



N 



(yi) 



(A3) 



Note that, in addition to the usual SP(2N) spin 
operators^ we have the density operator which is asym- 
metric in the fermion's species index: 

r b = Urjf a - /M, s« b3 = l(/ T Qt / T b - /f/f), 



gab+ — _ f fd\ t b | fb\ 



(f?fi + f?f!), s ab - = {s ab+ y 



(A4) 



The Hilbert space of this generalized spin Hamiltonian 
is defined as the SU (2)-gauge invariant subspace of the 
fcrmionic Hilbert space* 6 " 9 - That is, any fermion wavefunc- 
tion which respects the following local constraints is an 



[S ab3 ,S cd3 ] 

[s ab3 ,s cd+ ] 

[s ab3 ,s cd -] 
[s ab+ ,s cd+ ] 
[s ab+ ,s cd -] 

[i) ab ,S cd± ] 

[?p ab , s cd3 ] 



_^bcgad3 _ ^ad gcb3~^ 

= l(s bc S ad+ + 5 bd S ac+ ), 

= -l(s ad s bc - + s ac s bd -), 
= o, [s ab -,s cd -] = 0, 

^ jac <^bd3 _|_ ^ad g^ c 3 _|_ g&dZ _|_ ^bd gac3 j 

_ QQ-d^ 1 J ac ^j^dEi ^&d gcb± _|_ ^bd <£j ca 

^J^ c gadZ _ jac gjbdS ^ad g c ^3 _|_ fibd Q ca ^^j 

= l(5 bc tp ad - 6 ac ip bd ~ 5 ad ^ cb + S bd ^ ca ). 

(A5) 



Using them, one can argue that the generalized spin 
Hamiltonian given in Eq. (|A3|) is invariant under those 
continuous symmetries which are generated by 



N 



N 



N 



/ ab \ ^ naa3 \ y naal \ ^ naa2 
Ytot ' / j °tot i / j J tot ' / j °tot • 



a=l 



a=l 



When N = 1, V' ab disappears by itself and Eq. (|A3|) in 
combination with Eq. (|A5|) reduces to the SU (2) Heisen- 
berg spin model defined in Eq. |2]). 
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